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Abstract 

The Z„ elliptic Gaudin model with integrable boundaries specified by generic non- 
diagonal K-matrices with n + 1 free boundary parameters is studied. The commuting 
families of Gaudin operators are diagonalized by the algebraic Bethe ansatz method. 
The eigenvalues and the corresponding Bethe ansatz equations are obtained. 
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1 Introduction 



In the study of one-dimensional many-body systems with long-range interactions Gaudin 
type models [1] stand out as a particularly important class due to their applications in many 
branches of physics such as the BCS theory [2] of small metaUic grains [3, 4, 5], theoretical 
nuclear physics [6, 7], quantum chromo-dynamics (QCD) theory [8, 9] and Seiberg-Witten 
theory of supersymmetric gauge theory [10]. They also provide a powerful way for the 
construction of integral representations of solutions to the Knizhnik-Zamolodchikov (KZ) 
equation [11, 12, 13, 14, 15]. 

Rational and trigonometric Gaudin magnets have been extensively investigated in the 
literature. It is well-known that the periodic Gaudin's magnet Hamiltonians (or Gaudin 
operators) can be constructed via the quasi-classical expansion of the transfer matrix {row- 
to-row transfer matrix) of an inhomogeneous spin chain with periodic boundary conditions 
[16, 17]. Gaudin models with non-trivial boundary conditions can be in principle treated 
by means of the Sklyanin's boundary inverse scattering method [18]. In this case the quasi- 
classical expansion of the corresponding boundary transfer matrix (double-row transfer ma- 
trix) produces generalized Gaudin Hamiltonians with boundaries specified by certain K- 
matrices [14, 19]. In particular, twisted boundary conditions and open boundary conditions 
associated with diagonal K-matrices give rise to Gaudin magnets in non-uniform local mag- 
netic fields [14] and interacting electron pairs with certain non-uniform long-range coupling 
strengths [4, 20, 21, 22, 19]. Very recently, the XXZ Gaudin model with generic integrable 
boundaries specified by generic non-diagonal K-matrices [23, 24] was solved by the algebraic 
Bethe ansatz method [25]. 

So far, study of Gaudin type models via the algebraic Bethe ansatz method is quite 
advanced for rational and trigonometric interaction cases. On the other hand, the elliptic 
interaction case is rather less developed. We are aware of only two works [17, 15] where the 
XYZ Gaudin model and its face-type counterpart were respectively constructed and solved. 
In this paper we will solve, using the algebraic Bethe ansatz method, the most generic 
case-the Z„ elliptic Gaudin magnets with open boundary conditions specified by the generic 
non-diagonal K-matrices given in [26, 27], which contain n-\-l free boundary parameters. In 
section 3, we construct the elliptic Gaudin operators associated with the generic boundary 
K-matrices. The commutativity of these operators follows from the standard procedure [16, 
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14, 19] specializing to the inhomogeneous Z„ Belavin model with open boundaries [28, 29], 
thus ensuring the intcgrabihty of the Z„ elliptic Gaudin magnets. In section 4, we diagonalize 
the Gaudin operators simultaneously by means of the algebraic Bethe ansatz method. This 
constitutes the main new result in this paper. The diagonalization is achieved by means of 
the technique of the "vertex- face" transformation. In section 5, we conclude this paper with 
a summary and comments. 



2 Preliminaries: the inhomogeneous Belavin model 
with open boundaries 

Let us fix a positive integer n > 2, a complex number r such that Im{T) > and a generic 
complex number w. Introduce the following elliptic functions 
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{u,t) — exp {-\/^7r [(m + aj^T + 2(m + a) (w + 6)] } 
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{u, nr), a{u) — 9 



1 -\ 

2 
1 
2 



(2.1) 

(2.2) 
(2.3) 



Among them the cr-function^ satisfies the following identity: 

a{u + x)a{u — x)a{v + y)(j{v — y) — cr{u + y)(T{u — y)a{v + x)a{v — x) 
— a{u + v)a{u — v)a{x + y)(T{x — y). 

Let g, h, he n X n matrices with the elements 

hij = 5i+ij, Qij = co'-Sij, with CO = e~^, i,j e Z„. 

For any a = {ai,a2), cti, 0:2 G Z„, one can introduce an associated n x n matrix 1^ defined 

by 



and an elliptic function aa{u) given by 

aa{u) = 9 



1 I "1 

2 n 

1 I 02 

2 n -I 



(w,r), and o-(o,o)(m) = <y{u) 



^Our (T- function is the t?-function [30]. It has the following relation with the Weierstrassian a- 

function if denoted it by aw{u): aui{u) ex e^^"^a{u), rji = 4^^^ i-q'i^ ) andg = e^^"^ . Consequently, 

our ^-function given by (2.3) is different from the Weierstrassian C-function by an additional term — 2r/iu. 
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The Z„ Belavin R-matrix is given by [28] 

R-M = -^^^ J: '-^^^I. « (2.4) 

a(u + w) ^ naJ-) 

The R-matrix satisfies the quantum Yang-Baxter equation (QYBE) 

Rl2{ui - U2)Rl3{ui - U3)R23{U2 - U3) = i^23(^i2 " U3)Ri3{ui - U3)Ri2{ui - U2), (2.5) 

and the properties [31], 

Unitarity : Rf2{u)R^^{-u) = id, (2.6) 

Crossing-unitarity : (i?'')2\(-w - ntuMRy^Ju) = (7{u)cr{u + nw) ^ .^^ 

a[u + w)a{u + nw — w) 

Quasi-classical property : i?^(ii)|.u,_>o = id. (2.8) 

Here R^i{u) — -pL2-Ril(^i)^i2 with P12 being the usual permutation operator and ti denotes 
the transposition in the z-th space. Here and below we adopt the standard notation: for 
any matrix A e End(C"), Aj is an embedding operator in the tensor space C" (8) C" (8) • • •, 
which acts as A on the j-th space and as an identity on the other factor spaces; Rij{u) is an 
embedding operator of R-matrix in the tensor space, which acts as an identity on the factor 
spaces except for the i-th and j'-th ones. The quasi-classical properties (2.8) of the R-matrix 
enables one to introduce an associated Z„ classical r-matrix r{u) 



R^{u) = id^wr{u)^0{w'^), when w — ^ 0, 

_ ^"'^/-^ ^ , (T'(0)(7a(t() 1 , d 



r[u] 



E .'(0) = (2.9) 



In the above equation, the elliptic ^-function is defined in (2.3). 

One introduces the "row-to-row" monodromy matrix T{u) [33], which is an n x n matrix 
with elements being operators acting on (C")®^ 

T{u) = Rg{u + Zi)Rg{u + Z2)--- R^A'^ + zn). (2.10) 

Here {zi\i = 1, - ■ ■ , N} are arbitrary free complex parameters which are usually called inho- 
mogeneous parameters. With the help of the QYBE (2.5), one can show that T{u) satisfies 
the so-called "RLL" relation 

Rg{u - v)T,{u)T2{v) = T2{v)T,{u)Rg{u - v). (2.11) 
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An integrable open chain can be constructed as follows [18]. Let us introduce a pair of 

K- matrices K^{u) and K^{u). The former satisfies the refiection equation (RE) 

i?f2(lti - 1t2)i^r(ltl)i?fi(ltl + U2)K2{U2) 

= K^{u2)RUui + U2)K^iui)Rg{ui - U2), (2.12) 

and the latter satisfies the dual RE 

Ri2{u2 - Ui)K^{ui)R2i{-Ui -U2- nw)K^{u2) 

= K^{u2)R^2{-^i -U2- nw)K+{ui)R^^{u2 - Ui). (2.13) 

For the models with open boundaries, instead of the standard "row-to-row" monodromy 
matrix T{u) (2.10), one needs the "double-row" monodromy matrix T{u) 

T{u) = T{u)K-{u)T-\-u). (2.14) 

Using (2.11) and (2.12), one can prove that T{u) satisfies 

RuiUl - U2)Ti{ui)Rg{Ui + U2)T2{U2) = T2{u2)Ri2{ui + U2)Ti{ui) Rg{ui - U2). (2.15) 

Then the double-row transfer matrix of the inhomogeneous Z„ Belavin model with open 
boundary is given by 

r{u) = tr{K+{u)T{u)). (2.16) 

The commutativity of the transfer matrices 

[t{u),t{v)]^0, (2.17) 

follows as a consequence of (2. 5)- (2. 7) and (2.12)-(2.13). This ensures the integr ability of 
the inhomogeneous Z„ Belavin model with open boundary. 



3 Zn elliptic Gaudin models with generic boundaries 

In this paper, we will consider a non-diagonal K- matrix K~[u) which is a solution to the 
RE (2.12) associated with the Z„ Belavin R-matrix [26] 

^"(«)« = E l[m+u) (3-1) 

5 



The corresponding dual K- matrix K~^{u) which is a solution to the dual RE (2.13) has been 
obtained in [27]. With a particular choice of the free boundary parameters with respect to 
K~{u), we have 

""^"^^ - hill -(^-^^) N(A.+e— ^) 

x0i:U(-«)^i!Ui(^)- (3-2) 

In (3.1) and (3.2), 0, 0, are intertwiners which will be specified in section 4. We consider 
the generic {Aj} such that Aj 7^ Xj {modulo Z + rZ) for i 7^ j. This condition is necessary for 
the non-singularity of K^{u). It is convenient to introduce a vector A = Y17=i -^i^i associated 
with the boundary parameters {Aj}, where {e^, ^ = 1, • • • , n} is the orthonormal basis of the 
vector space C" such that (ej,ej) = Sij. 

As will be seen from the definitions of the intertwiners (4.8), (4.10) and (4.11) specialized 
to m = A, (t>\,x-wi{u) and 4>\,\-wiiu) do not depend on w but (i)\\-wi{u) does. Consequently, 
the K- matrix K''{u) does not depend on the crossing parameter to, but K^{u) does. So we 
use the convention: 

Kiu) = lim K- iu) = K- (u) . (3.3) 

We further restrict the complex parameters ^ and ^ to be the same, i.e. 

l=i, (3.4) 

so that (3.5) below is satisfied. Hence, the K-matrices depend on n + 1 free parameters 
= 1, • • • ,n} and ^, which specify the integrable boundary conditions [24]. Moreover, 
the K-matrices satisfy the following relation thanks to the restriction (3.4) 

\im{K+(u) K-{u)} = lim {K+{u)}K{u) = id. (3.5) 

Let us introduce the Z„ elliptic Gaudin operators {Hj\j = 1, 2, • • • , N} associated with the 
inhomogeneous Z„ Belavin model with open boundaries specified by the generic K-matrices 
(3.1) and (3.2): 

TV f ^ 1 

= r^(%) + Yl ^^i^^i ~ + { Yl ^i^^^i + f ^i(^j)' (3-6) 
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where Tj{u) = £{K,{u)}\.^=oKjiu), j = 1,---,7V, with Kj{u) = tro {K+{u)R^j{2u)Poj}. 
Here {zj} are the inhomogeneous parameters of the inhomogeneous Z„ Belavin model and 
r{u) is given by (2.9). For a generic choice of the boundary parameters {Ai, •••,A„, 
Tj{u) is a non-diagonal matrix. 

The elliptic Gaudin operators (3.6) are obtained by expanding the double-row transfer 
matrix (2.16) at the point u — zj around w — 0: 

T{zj) = T{zj)\^=o + wHj + 0{w^), j = l,---,N, (3.7) 
Hj = -^T{zj)U=o. (3.8) 

The relations (2.8) and (3.5) imply that the first term T{zj)\w=o in the expansion (3.7) is 
equal to an identity, namely, 

r(z,)|,»=o = id. (3.9) 

Then the commutativity of the transfer matrices {T{zj)} (2.17) for a generic w implies 

[Hj,Hk] = 0, i,j = l,---,N. (3.10) 

Thus the elliptic Gaudin system defined by (3.6) is integrable. Moreover, the fact that 
the Gaudin operators {Hj} (3.6) can be expressed in terms of the transfer matrix of the 
inhomogeneous Z„ Belavin model with open boundary enables us to exactly diagonalize 
the operators by the algebraic Bethe ansatz method with the help of the "vertex-face" 
correspondence technique, as will be shown in the next section. The aim of this paper is to 
diagonalize the elliptic Gaudin operators Hj, j — 1, ■ ■ ■ , N (3.6) simultaneously. 

4 Eigenvalues and Bethe ansatz equations 

4.1 A^^li SOS R- matrix and face- vertex correspondence 

The An-i simple roots {ai} can be expressed in terms of the orthonormal basis {ej} as: 

Q!j = ej-ej+i, i = 1, • • • ,n - 1, 
and the fundamental weights {Aj | i = l,---,n — 1} satisfying (Aj, Uj) — Sij are given by 

i . n 

Ai = 5^efc-^ J]efc. 



n 

k=l k=l 



Set 

^ n n 

e, e=— ^^ejfc, i — l,---,n, then — 0. (4-1) 



^ k=l i=l 



For each dominant weight A = ^^I^ cii^i, di £ Z+ (the set of non- negative integers), there 
exists an irreducible highest weight finite-dimensional representation Va of v4„_i with the 
highest vector |A). For example the fundamental vector representation is Va.i- 

Let 1^ be the Cartan subalgebra of A^-i and i)* be its dual. A finite dimensional diag- 
onalizable [)-module is a complex finite dimensional vector space W with a weight decom- 
position W = ®^^^*W[fi], so that [) acts on W[fi] by xv = fi{x)v, {x E i), v E W[ij]). 
For example, the fundamental vector representation Vaj = C", the non- zero weight spaces 
W[i] =Cei, i = l,---,n. 

For a generic m e C", define 

rrii = (m, e^), rriij = - rnj = {m, - e^), i,j = l,---,n. (4.2) 
Let R{u, m) e EndiU" ® C") be the R-matrix of the ^J^i SOS model [34] given by 

n 

R(u, m)^Yl 0Eii + J2 rn)Eii Ejj + E^jiu, m)Eji ® Eij} , (4.3) 

i=l iytj 

where E^j is the matrix with elements {Eij)\ — SjkSu. The coefficient functions are 

KM^l. (4.4) 

Ri;^,m) = ^P^^^^, (4.6) 
■' a[u + w)a[mij) 

and rriij are defined in (4.2). The R-matrix satisfies the dynamical (modified) QYBE 

Ri2{ui -U2,m- wh^^^)Ris(ui - u^,m)R2z{u2 -Us,m- wh^^^) 

= R2z{u2 - U3,m)Ris{ui -U3,m- wh^'^^)Ri2{ui - U2,m), (4.6) 

and the quasi-classical property 

R{u,m)\^^o = id. (4.7) 

We adopt the notation: Ri2{u,m — wh^^^) acts on a tensor vi i>2 <8) as R{u,m — 'Wfj,)<^idif 
V3 e Moreover, the R-matrix satisfies the unitarity and the modified crossing-unitarity 

relation [34]. 



Let us introduce an intertwiner — an n-component column vector 0m,m-iDi(w) whose k-th 
element is 

€,„ji^)-0^'\^ + nm,). (4.8) 

Using the intertwining vector, one derives the following face- vertex correspondence relation 
[34] 

= J2 Riui - U2, m)fj (l)„i^„^i+k){ui) (8) 0^,„K«2). (4.9) 

k,l 

Then the QYBE of Z„ Belavin's R-matrix R^(u) (2.5) is equivalent to the dynamical Yang- 
Baxter equation of A^^li SOS R-matrix R{u,m) (4.6). For a generic m, we may introduce 
other types of intertwiners (f), (f) satisfying the conditions, 

n 

E Cm-.o{u) = V, (4.10) 

fe=l 
n 

fc=l 



from which one derives the relations. 



n 



E <t^^-n,iSy) = (4.12) 



n 



E'^2+.A,n.W'^!^Un.N = ^^.-- (4-13) 

With the help of (4. 10)- (4. 13), we obtain the following relations from the face- vertex corre- 
spondence relation (4.9): 



m+wk,m 



(Ui) (g) idj -Rf2(til - U2) (id (g) (f>m+wj,m{u2)^ 
= E ^(^1 ~ m)g 0m+«;(J+i),m+«;j(Ml) ® 0^+^(fc+f),^+^fe(M2) , (4.14) 



i,l 



(j>m+wk,miUi) ® '^m+«;(fe+0,Wfe(^2)) Rui^l " ^2) 

= E ^(""1 ~ '"2, "^)»^i 0m+«;(i+i),m+«;i(Ml) 0m+«;i,m('"2) , (4.15) 



id ® 0m,m-«;[("2)) ^f2(lil " ^2) (^0m,m-m(Ml) O id^ 



= R{ui - U2, m)g (p„^-^l,m-u,(k+i)i'^l) ® (Pm-wi,m-w{i+j){u2) , (4.16) 
= ^ R{Ui - U2, m)fj 4>m,m-wi{Ul) ® 4>m-wi,m-w{i+3){u2) ■ (4.17) 

The face-vertex correspondence relations (4.9) and (4.14)-(4.17) will play an important role 
in translating formulas in the "vertex form" into those in the "face form" so that the algebraic 
Bethe ansatz method can be applied to diagonalize the transfer matrix. 

Corresponding to the vertex type K-matrices (3.1) and (3.2), one introduces the following 
face type K-matrices /C and K,, as in [27] 

IC{X\u)i = ^l-^(.-i), A-..-(«) K{u)t ,_,(-«), (4.18) 

s,t 

IC{X\u)l = E 4:\-.i(-«) ^{u)t ,_,(«). (4.19) 

s,t 

Through straightforward calculations, we find the face type K-matrices simultaneously have 
diagonal forms^ 

}CiX\u)i = 5i kiX\u; Oi, iC{X\u)i = ~kiX\u)i, (4.20) 
where functions k{X\u;^)i, k{X\u)i are given by 

MA|n;a. = (4.21) 

Moreover, one can check that the matrices }C{X\u) and IC{X\u) satisfy the SOS type reflection 
equation and its dual, respectively [27]. Although the K-matrices K'^i^u) given by (3.1) and 
(3.2) are generally non-diagonal (in the vertex form), after the face- vertex transformations 
(4.18) and (4.19), the face type counterparts ]C{X\u) and ]C{X\u) simultaneously become 
diagonal. This fact enables the authors to apply the generahzed algebraic Bethe ansatz 
method developed in [29] for SOS type integrable models to diagonalize the transfer matrix 
t{u) (2.16). 

^As will be seen below (4.52), the spectral parameter u and the boundary parameter ^ of the reduced 

double-row monodromy matrices constructed from IC{X\u) will be shifted in each step of the nested Bethe 
ansatz procedure [29]. Therefore, it is convenient to specify the dependence on the boundary parameter ^ 
of )C{X\u) in addition to the spectral parameter u. 
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4.2 Algebraic Bethe ansatz 

By means of (4.12), (4.13), (4.19) and (4.20), the transfer matrix t{u) (2.16) can be recast 
into the following face type form: 

t{u) = tr{K+{u)T{u)) 

= Y,mu)>tr{\\ur^ = Y,~mu),r{\\u)>;. (4.23) 

Here we have introduced the face-type double-row monodromy matrix T(A|k), 

= <Px-w{il-i>),X-wti{u) T(?7,) (t)X,X-wij.{-u) 

This face-type double-row monodromy matrix can be expressed in terms of the face type 
R-matrix R{u, A) (4.3) and the K-matrix K:{X\u) (4.18) [29]. Moreover from (2.15), (4.9) and 
(4.13) one may derive the following exchange relations among T{X\u)1^: 



>12 



X R{u^ + xy^ r(A + w{h + h)\u2)t 

X r(A + «;(i2 + ^2)|«i)i^ R{u, - U2, Xf^. (4.25) 
As in [29], we introduce a standard notation for convenience^: 

A{\\u) = T{\\u)\, B,i\\u)^'^^^^, t = 2,...,n, (4.26) 

Vi{X\u) = - ^v^) ^^x\u)i diRi2u, A + wi){] A{X\u)} , 

^,J = 2, (4.27) 



^The scalar factors in the definitions of the operators B{X\u) and r'(A|u) are to make the relevant com- 
mutation relations as concise as (4.28)-(4.30). 
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Prom (4.25) one may derive the commutation relations among A{X\u), V{X\u) and B{X\u). 
Here we give those which are relevant for our purpose 



Bi{X + w{i - i)\v)AiX + w{i - i)\u) 

Bi{X + w{i-i)\u)A{X + w{i-i)\v) 



A{X\u)Bi{X + w{i-l)\v) 

a{u + v)a{u — V — w) 
a{u + V + w)a{u — v) 

a{w)a{2v) a{u — v — Xii + w) 
a{u — v)a{2v + w) cr{Xii — w) 

a(w) yr^ criu + V + Xai + 2w) ^ ,^ ,^ , 
7 — T 7 ^ 7^ —^ -Ba{X + w{a-l) u)V^{X + w{i-l) v), 

^ ' 06=2 ^ ' 

(4.28) 

Vf{X\u)B,iX + wij-i)\v) 

a{u — V + w)a{u + v + 2w) 
a{u — v)a{u + V + w) 

X I j2 Riu + v + w,X- wit^Riu -v,X + wjfr 

Kai,a2,l3i, 132=2 

X B^,{X + w{k + h-i- l)\v)Vll{X + w{3 - l)\u)] 
a{w)a{2u + 2w) J a{u — v + Xia — w) 



a{u- v)a{2u + w) ) ^ a{Xia - w 



X Bf^{X + w{k + p-i- i)\u)Vf{X + w{j - i)\v) } 



_^ a{w)a{2v)a{2u + 2w) 



X 



a{u + V + w)a{2v + w)a{2u + w) 

R(2u + w,X- m)^f 

X Ba{X + w{k + a-i- i)\u)A{X + wQ - i)\v)^ , (4.29) 
Bi{X + w{i - i) \u)Bj{X + w{i + j- 2i)\v) 

n 

= J] i?(M-'y,A-2wi)ff S^(A + 'u;(/3- i)|^;)S„(A + 'u;(a + /3-2i)|M). (4.30) 

a,P=2 

In order to apply the algebraic Bethe ansatz method, in addition to the relevant com- 
mutation relations (4.28)-(4.30), one needs to construct a pseudo- vacuum state (also called 
reference state) which is the common eigenstate of the operators A, V\ and is annihilated by 
the operators Cj. In contrast to the trigonometric and rational cases with diagonal K^{u) 
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[18, 35], the usual highest- weight state 



/ 1 \ 





/1\ 





V ^ / \ - J 

is no longer the pseudo-vacuum state. However, after the face- vertex transformations (4.18) 
and (4.19), the face type K-matrices )C{\\u) and IC{\\u) simultaneously become diagonal. 
This suggests that one can translate the Z,„ Belavin model with non-diagonal K-matrices 
into the corresponding SOS model with diagonal K-matrices IC{X\u) and IC{\\u) given by 
(4.18)-(4.19). Then one can construct the pseudo- vacuum in the "face language" and use 
the algebraic Bethe ansatz method to diagonalize the transfer matrix [29]. 
Let us introduce the corresponding pseudo- vacuum state \ fl) 



■'X-{N-l)wi,X-Nwi^ 



^A-(iV- 



-2)«,i,A-(iV-l)«;i(-^2) • • • (8 </'A,A-«;i(-^iv). (4.31) 



This state depends only on the boundary parameters {Aj} and the inhomogeneous parameters 
{zj}, not on the boundary parameter C^. We find that the pseudo- vacuum state given by (4.31) 
satisfies the following equations, as required. 



AiX - Nwi\u)\Q) 

vi{x - Nwi\u)\n) 



kiX\u■0^\Q), 

5i/W(«)fc(A|« + |;e-|), 



X ■ 



N 

n 



a{u + Zk)a{u - Zk) 



k=i + + w)a(u — Zk + w) 
Ci{X - Nwi\u)\Q,) = 0, z = 2,---,n, 
Bi{X - Nwi\u)\n) ^ 0, i = 2, 



,n. 



(4.32) 



^,J =2,- • -,71,(4.33) 

(4.34) 
(4.35) 



Here /'•^•'(■u) is given by 



a{2u)a{Xi + u + w + 



(4.36) 



a{2u + w)a{Xi + u + C)' 

In order to apply the algebraic Bethe ansatz method to diagonalize the transfer matrix, 
we need to assume N — n x I with I being a positive integer [29]. For convenience, let us 
introduce a set of integers: 



Ni ^ {n - i) X I, i = 0, 1, • • • , n - 1, 



(4.37) 
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and complex parameters {v^^''\ k = 1,2, ■ ■ ■ , A^j+i, i = 0,1, ■ ■ ■ ,n — 2}. As in the usual 

nested Bethe ansatz method [36, 37, 35, 38, 29], the parameters {v^^^} will be used to specify 
the eigenvectors of the corresponding reduced transfer matrices. They will be constrained 
later by the Bethe ansatz equations. For convenience, we adopt the following convention: 

Vk = v'^^\ k=l,2,---,N,. (4.38) 

We will seek the common eigenvectors (i.e. the so-called Bethe states) of the transfer matrix 
in the form 

n 

\vi,---,VN,) = ^''''''"■■''^i'Bi,(A + «;(?i-i)|t;i)Bi,(A + «;(£i + Z2-2i)|^;2) 

n,---i«jvi=2 

Ni 

X • • • Bi^^ {X + wJ2ik- wNii\vN, m. (4.39) 

k=l 

The indices in the above equation should satisfy the following condition: the number of 

ik = j, denoted by #(j), is /, i.e. 

#(j) = /, j = 2,...,n. (4.40) 

With the help of (4.23), (4.26) and (4.27) we rewrite the transfer matrix (2.16) in terms 
of the operators A and T>1 

n 

r{u) = J2~^{\\u),r{\\u): 

n 

= fc(A|ti)M(A|ti) + ^fc(A|ii),r(A|ii)^ 

i=2 
n 

= k{\\uy A{\\u) + ^ k{\\u)i R(2u, A + wi){] A{X\u) 

i=2 
n 

+ k{\\u)i (r{\\u)\ - R{2u, A + wi)\\A{X\i 

i=2 

n 

= YHMu)iR{2u,X + wi){\A{\\u) 
1=1 

+ g fc(^) {X\u + |), ^^^^^ (T(A|«)i - R{2u, A + wl)\\A{X\u)) 

n 

= a^'\u)A{\\u) + J2k^'\X\u+^),V{Muyi. (4.41) 

i=2 
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Here we have used (4.27) and introduced the function a'^^^u), 

n 

= ~k{\\u)i R(2u, X + wi)\], (4.42) 

i=l 

and the reduced K-matrix K,^^\X\u) with the elements given by 

}C^'\X\u)l = dlk^'HXlu}^, ^,j = 2,---,n, (4.43) 

To carry out the nested Bethe ansatz process [36, 37, 35, 38, 29] for the Z„ Belavin model with 
the generic open boundary conditions, one needs to introduce a set of reduced K-matrices 
}C^"^\X\u) [29] which include the original one IC{X\u) = }C^^\X\u) and the ones in (4.43) and 
(4.44): 

^M(A|m)^' = 6l¥"'\X\u)^, i,j = m + l,---,n, m = 0,---,n-l, (4.45) 

Moreover we introduce a set of functions {a^'^\u)\m — 1, • • • ,n — 1} (including the one in 
(4.42)) related to the reduced K-matrices K,^'^\X\u) 

n 

Q;M(ii) = ^i?(2u,A + «;m)j;^^("^-^)(A|u)i, m = l,---,n. (4.47) 

i=m 

Carrying out the nested Bethe ansatz, we find [29] that with the coefficients F'i>'2> - >«Afi jn 
(4.39) properly chosen, the Bethe state \vi, • • • , v^^) is the eigenstate of the transfer matrix 
(2.16), 

t{u)\vi,---,vn,) = A{u; ^, {vk})\vi, ■ ■ ■ ,vni) , (4.48) 
with eigenvalues given by 

cr{Xi + ^ + u) a{u + Vk + w)a{u - Vk) 
a{2u)a{Xi -\- u-\-w + ^) ^ a{u — Vk-\- w)a{u -\-Vk-\- 2w) 

rr(9.ii -\- nArr(\-. -\- ii -I- H 



a{2u + w)a{Xi + + ~ Vk)cr{u + Vk + w) 

No 

X 

k 



ft , + Aa)(. + ^; ^ - ^, {.W}).(4.49) 

-'-i(T(M + ^fe + 'u;)(T(ii-Zfe + «;) ^ 2'^ 2'^ ^ ^ 
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The eigenvalues {A^'^\u; {v^k^})} (with the original one A{u; {vk}) = A^^^\u; ^, {v^k^})) of 
the reduced transfer matrices are given by the following recurrence relations 



a{2u)(T{X,+i + u + W + ^(^)) yr a{u - + w)a{u + + 2w) 

X n + f -f . 

,^^ia{u + zl' +w)a{u- zl' +w) ^ ^ 

i = l,---,n-2. (4.50) 

The reduced boundary parameters {^'•*-'} and inhomogeneous parameters {-2^*''} are given by 

e^-^^ = e« - f , 4'"^^ = + f , i = 0, . . . , n - 2. (4.52) 

Here we have adopted the convention: ^ = 4°'' = -^fe- The complex parameters {4*"*} 
satisfy the following Bethe ansatz equations: 

a{Xi + C - Vs)a{2vs + w) ^ a{vs + Vk)a{vs - Vk - w) 



(j(Ai + C + + w)(t(2i;s + 2w) ^ (^(v^ + f + 2w)a{vs - Vk + w) 



N 



TT , ^(^- + ^^)^(^--^^) ^ A(i)(,^ + ^; e - {4^^}), (4.53) 
l}^a{vs + Zk + w)a{vs-Zk + w) ^ « 2 ' 2 ' ^ ^ ^ ' 

a(A,+i + - v^^)ai2vi'^ + w) a(4') + 4'V(4^^ - 4'^ - w) 



' a(A,+i + C« + 4'^ + w)a{2vf^ + 2w),4,Li ^(4'^ + vf + 2w)(t(4'' - vf + u;) 

TT a(t;. )a(^;. - ) ^(i+i). « _ ^ r (i+i)-. ^ 

11 / (i) , (j) , X / (i) (j) , X V^s ^ 9>'S 9>\^fe 



i = l,---,n-2. (4.54) 

4.3 Eigenvalues and associated Bethe ansatz equations 

The relation (3.8) between {Hj} and {T{zj)} and the fact that the first term on the r.h.s. of 
(3.7) is a c-number enable us to extract the eigenstates of the elliptic Gaudin operators and 
the corresponding eigenvalues from the results obtained in the previous subsection. 
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Introduce the functions {f3^'^\u,w)} 

p(^-^^\u,w)^P'''-'\u)^a^^^^\u) "^^'^'^i'T'""^ V i-0,.--,n-2. (4.55) 

o-(Aj+i + t, + u + w - |w) 

Then by (4.50), (4.51) and (4.52), the Bethe ansatz equations (4.53) and (4.54) become, 
respectively, 

^ ^ ' ^2^M-(-in.mf)a(.«-.mf) 

i = 0,---,n-3, (4.56) 

a(2^r^) + w) <7(^;r'^ + v^Mv^-'^ - - w) 



X 



Here we have used the convention: Vf, — Zk, k — 1, ■ ■ ■ , N . The quasi-classical property 
(4.7) of R{u,m), (4.47) and (4.55) lead to the following relations 

P^'+^\u,0)^l, ■^(3^'+^\u,0) ^0, i^0,---,n-2. (4.58) 

ou 

Then, one may introduce functions {7*^*"'"^)(ti)} associated with {P^^^^\u,w)} 

l^'-''\u) = -^P^'^'\u,w)\^=0: i = 0,.-.,n-2. (4.59) 

Using (4.58), wc can extract the eigenvalues hj (resp. the corresponding Bethe ansatz equa- 
tions) of the Gaudin operators Hj (3.6) from the expansion around w = for the first order 
of w of the eigenvalues (4.49) of the transfer matrix t{u — Zj) (resp. the Bethe ansatz 
equations (4.56) and (4.57) ). Finally, the eigenvalues of the Z„ elliptic Gaudin operators 
are 

Ni 

hj = l^'\zj) + C(Ai + C + Zj)-Yl + ^k) + C{zj - Xk)} , (4.60) 

k=l 
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where C-function is defined in (2.3). The parameters k = 1,2, ■ ■ ■ , Ni+i, i = 

0, 1, ■ ■ ■ , n — 2} (inchiding Xk as Xk = xf \ k = 1, • ■ • , Ni) are determined by the following 
associated Bethe ansatz equations 



7^ \x 



s 



k^s,k=l 



k=l 

Ni+2 



k=l 

i = Q,---,n-3, (4.61) 



JV„_1 

k^s,k=l 



AT, 



n-2 



- E {c(4"-^^ + 4""'^) + cix^r'^ - xf^)} . (4.62) 
fe=i 

Here we have used the convention: x\^^^ — Zk, k — 1, ■ ■ ■ , N in (4.61). 

5 Conclusions 

We have studied the Z„ elliptic Gaudin model with generic boundaries specified by the non- 
diagonal K- matrices K^{u), (3.1) and (3.2). In addition to the inhomogeneous parameters 
{zj}, the associated Gaudin operators {Hj}, (3.6), have n + 1 free parameters {A^} and 
^. As seen from section 4, although the "vertex type" K-matrices K^{u) (3.1) and (3.2) 
are non-diagonal, the compositions, (4.18) and (4.19), lead to the diagonal "face-type" 
K-matrices after the face-vertex transformation. This enables us to successfully construct 
the corresponding pseudo-vacuum state \Q) (4.31) and apply the algebraic Bethe ansatz 
method to diagonalize the transfer matrix t{u) of the inhomogeneous Z„ Belavin model 
with generic open boundaries. Furthermore, we have exactly diagonalized the generalized 
Gaudin operators {Hj}, and derived their eigenvalues (4.60) as well as the associated Bethe 
ansatz equations (4.61) and (4.62). Taking the scahng hmit [39] of our general results for 
the special n — 2 case, we recover the results obtained in [25] for the XXZ Gaudin model 
with generic open boundaries. 
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